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Abstract 

A chain of metallic particles, of sufficiently small diameter and spacing, allows linearly polarized 
plasmonic waves to propagate along the chain. In this paper, we consider how these waves are 
altered by an anisotropic host (such as a nematic liquid crystal) or an applied magnetic field. In a 
liquid crystalline host, with principal axis (director) oriented either parallel or perpendicular to the 
chain, we find that the dispersion relations of both the longitudinal (L) and transverse (T) modes 
are significantly altered relative to those of an isotropic host. Furthermore, when the director 
is perpendicular to the chain, the doubly degenerate T branch is split by the anisotropy of the 
host material. With an applied magnetic field B parallel to the chain, the propagating transverse 
modes are circularly polarized, and the left and right circularly polarized branches have slightly 
different dispersion relations. As a result, if a linearly polarized transverse wave is launched along 
the chain, it undergoes Faraday rotation. For parameters approximating that of a typical metal 
and for a field of 2T, the Faraday rotation is of order 1° per ten interparticle spacings, even taking 
into account single-particle damping. If B is perpendicular to the chain, one of the T branches 
mixes with the L branch to form two elliptically polarized branches. Our calculations include 
single-particle damping and can, in principle, be generalized to include radiation damping. The 
present work suggests that the dispersion relations of plasmonic waves on chain of nanoparticles 
can be controlled by immersing the chain in a nematic liquid crystal and varying the director axis, 
or by applying a magnetic field. 

PACS numbers: 78.67.Bf,64.70.pp,78.20.Ls,73.20.Mf 
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I. INTRODUCTION 



The optical properties of small metal particles have been of interest to physicists since 
the time of Maxwell|l|. Such particles, if subjected to light of wavelength much larger than 
their linear dimensions, exhibit optical resonances due to localized electronic excitations 
known as "particle" or " surface" plasmons. These plasmons can give rise to characteristic 
absorption peaks, typically in the near-infrared or the visible, which may play an important 
role in the optical response of dilute suspensions of metal particles in a dielectric hostj2-4|. 

Because of recent advances in sample preparation, it has become possible to study ordered 
arrays of metal particles in a dielectric host|5|-l8|. In one- dimensional ordered arrays of such 
closely-spaced particles, waves of plasmonic excitations can propagate along the chains. 



provided that the interparticle spacing is small compared to the wavelength of light 



20|. In this limit, the electric fields produced by the dipole moment of one nanoparticle 



induces dipole moments on the neighboring nanoparticles. The dispersion relations for both 
transverse (T) and longitudin al i V) plasmonic waves can then be calculated in the so-called 



quasistatic approximation 



10l-|l2l|. in which the curl of the electric field is neglected. While 



this approximation neglects some significant coupling between the plasmonic waves and free 
photons jiol. it gives reasonable results over most of the Brillouin zone. Interest in such 
plasmonic waves has grown vastly in recent years 12, isl . 

In this paper, we extend the study of propagating plasmonic waves in two ways. First, 
we calculate the dispersion relation for such plasmonic waves when the metallic chain is 
immersed in an anisotropic host, such as a nematic liquid crystal (NLC). Using a simple 
approximation, we show that both the L and T waves have modified dispersion relations 
when the director is parallel to the chain axis. If the director is perpendicular to that axis, 
we show that the previously degenerate T branches are split into two separate branches. 
Secondly, we consider the effects of a static magnetic field applied either parallel and per- 
pendicular to the chain. For the parallel case we show that a linearly polarized T wave is 
rotated as it propagates along the chain. For a field of 2 tesla and reasonable parameters for 
the metal, this Faraday rotation may be as large as 1-2° over ten interparticle spacings. A 
perpendicular field mixes together the L branch and one of the T branches, leading to two 
elliptically polarized branches. These results suggest that either an NLC host or an applied 
magnetic field could be used as an additional "control knob" to manipulate the properties 
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of the propagating waves in some desired way. 

The remainder of this paper is organized as follows. In the next section, we present the 
formalism which allows one to calculate the dispersion relations for L and T waves in the 
presence of either an anisotropic host or an applied dc magnetic field. In Section III, we give 
simple numerical examples, and we follow this by a brief concluding discussion in Section 
IV. 



II. FORMALISM 



A. Overview 



We consider a chain of identical metal nanoparticles, each a sphere of radius a, arranged 
in a one-dimensional periodic lattice along the z axis. The n*^ particle is assumed centered 
at (0,0, rac?) (—00 < n < +00). The propagation of plasmonic waves along such a chain 
of nanoparticles has already been considered extensively for the case of isotropic metal 



particles embedded in a homogeneous, isotropic medium |10l|. Various works have considered 
the quasistatic case in which the electric field is assumed to be curl-free; this is roughly 
applicable when both the radius of the particles and the distance between them are small 



compared the wavelength of light 
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121 ]. The extension of such studies to include radiative 
corrections, i. e., to the case when the electric fields cannot be approximated as curl-free, has 
also been carried out; these corrections can be very important even in some long- wavelength 
regimes 16 1. 



Here we consider how the plasmon dispersion relations are modified when the particle 
chain is immersed in an anisotropic dielectric, such as an NLC, or subjected to an apphed 
dc magnetic field. For the case of metallic particles immersed in an NLC, we assume that 
the host medium is a uniaxial dielectric, with principal dielectric constants e_L, e_L, and ey. 
For metal particles in the presence of an applied magnetic field, we take the host medium 
to be vacuum, with dielectric constant unity. 

In the absence of a magnetic field, the medium inside the particles is assumed to have 
a scalar dielectric function e{uj). If there is a magnetic field along the z axis, the dielectric 
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function of the particles becomes a tensor, whose components may be written 



e^y{uj) = -ey^(w) = iA{uj), (1) 



with all other components vanishing 



2l| . In the calculations below, we will assume that 



the nanoparticles are adequately described by a Drude dielectric function. In this case, the 



components of the dielectric tensor take the form 21 1 



= 1 - , ^ (2) 



and 

w (1 — iur) 

Here Up is the plasma frequency, r is a relaxation time, and Uc = eB/{mc) is the cyclotron 
frequency, where B = Bz is the magnetic field, m is the electron mass, and e is its charge. 
We will use Gaussian units throughout. In the limit ur — t- oo, we may write 

,2 

9 ' 



A{^) = #. (4) 



B. Uniaxially Anisotropic Host 

We first assume that the host has a dielectric tensor th with principal components e±, 
e±, and ey. Such a form is appropriate, for example, in a nematic liquid crystal below its 
nematic-to-isotropic transition. We begin by writing down the electric field at x due to a 
sphere with a polarization P(x'). In component form, this field takes the form (see, e. g., 



Ref 



22|) 



^.(x) = - j Gj^i^ - x')P,(x')rfV, (5) 

where repeated indices are summed over, and we use the fact that Qji = Qij. In eq. (I5l), 
P(x') = (e — e/i)E(x') is the polarization of the metallic particle, e is the dielectric function 
of the metal particle, and Q denotes a 3x3 matrix whose elements are 

^.,=9ia,G(x,-x;), (6) 
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where G(x — x') is a Green's function which satisfies the differential equation (see, e. g., Ref. 



22 



23|) 



V ■ e,, VG(x - x') = -(5(x - x') . (7) 

If the host dielectric tensor is diagonal and uniaxial with diagonal components e±, and 

ey, which we take for the moment to be parallel to the x, y, and z axes respectively, this 
Green's function is given by 22 1 



^(x 



1 



47re I e 



1/2 



{x - x'f + {y - y'Y , ( 



z — z 



2^ -1/2 



(8) 



Physically, — ^jj(x — x') represents the i component of electric field at x due to a unit point 
dipole oriented in the direction at x', in the presence of the anisotropic host. 

The next step is to use this result to obtain a self-consistent equation for plasmonic waves 
along a chain immersed in an anisotropic host. To do this, we consider the polarization of the 
n*'' particle, which we write as Pn(x) = 5eEj„ „(x), where Ej„„(x) is the electric field within 
the v^^ particle and 5e = e—eh- This field, in turn, is related to the external field acting on the 
n*'^ particle and arising from the dipole moments of all the other particles. We approximate 
this external field as uniform over the volume of the particle, and denote it 'Eiext,n- This 
approximation should be reasonable if the particle radius is not too large compared to the 
separation between particles (in practice, an adequate condition is probably a/d < 1/3, 
where a is the particle radius and d the nearest neighbor separation). Then Ej„ „ and 'E'ext,n 



are related by 



23| 



E,; 



ext,ni 



(9) 



where F is a "depolarization matrix" defined, for example, in Ref. 22|]. Eext,n is the field 
acting on the n*^ particle due to the dipoles produced by all the other particles, as given by 
eq. ([5]). Hence, the dipole moment of the n*'^ particle may be written 



Pr. 



a^tE 



ext,n 



where 



t = 6e(l-T6e 



-1 



(10) 



(11) 



is a "t-matrix" describing the scattering properties of the metallic sphere in the surrounding 
material. Finally, we make the assumption that the portion of Eea;f_„ which comes from 
particle n' is obtained from eq. ([5]) as if the spherical particle n' were a point particle located 



at the center of the sphere (this approximation should again be reasonable if a/d < 1/3). 
With this approximation, and combining eqs. ([5]), f lTOj) . and f lTTj) . we obtain the following 
self-consistent equation for coupled dipole moments: 



Pn 



-t ^ G{y^n - X„')Pn'. 



(12) 



Let us first assume that the principal axis of the anisotropic host coincides with the chain 
direction, which we take as the z axis. In this case, the L and T waves decouple and can be 
treated independently, because one of the principal axes of the T tensor coincides with the 
chain axis. First, we consider the longitudinally polarized waves. To find their dispersion 
relation, we need to calculate Qzz{^n — x„'). From the definition of this quantity, and from 
the fact that x„ = ndz = ZnZ, we can readily show that Qzzi^n ~ x„/) = ^ ^ 
Hence, we obtain the following equations for the Pnz's: 



Sen 



Pn 



3ex 



rii5ei 



Pn'. 



yl 13 ■ 



(13) 
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For transverse modes, the relevant Green's function takes the form Qxx{^n — ^n'] 
J iTTp"' '^^^ resulting equation for the dipole moments takes the form 



Pr. 



1 ^11 3 

3e^ 1 



5ei 



P-n 



Zr,'\ 



(14) 



rx5e± , , 

In the isotropic case with a vacuum host, ey = e_L = 1, and F^..^ = Tyy = Tzz = —1/3, The 
equations for both the parallel and perpendicular cases reduce to the results obtained in 
Ref. [10|] for both L and T modes, as expected. 

For an anisotropic host and only nearest neighbor interactions, the dispersion relation for 
the T waves is implicitly given by 



1 = - 



2a^ 



5e± 



3 rf3 1 - F I 5e I e2 



cos kd 



and for the L waves by 



6e 



1 



3rf3l-F||(5e|| ex 



cos kd. 



(15) 



(16) 



=11 I -| 

The forms of F^ and Fy are well known (see, e. g., Ref. |i22i], where they are denoted F^ 



and F-,-,). We rewrite them here for convenience: 



1 

1 
2 



sm 



^1 



Til + 



sm 



(17) 
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where \ = 1 — e±/e\\. If we assume that the metalhc particle has a Drude dielectric function 
of the form e{uj) = 1 — cjp/o;^, then the dispersion relation for T waves becomes 

T{k) 



u:l {l-e^)T{k)-V 



(18) 



where 



and that of the L waves is 



where 



T(fc)=r^-^^^cos(fcrf), (19) 



:i-e||)L(A;)-l' 



L(A;) = r|| + ^^— cos(fcrf). (21) 

o (X €_|_ 

Eqs. f|T8l) and fl20l) neglect damping of the waves due to dissipation within the metallic 
particles. To include this effect, one can simply solve eq. ( ITSl) or ( IT6l) for A;(ci;), using the 
Drude function with a finite r. The resulting fc(a;) will be complex in both cases; the inverse 
of the imaginary part of k{uj) will give the exponential decay length of the T or L wave 
along the chain. 

Now, let us repeat this calculation but with the principal axis of the liquid crystalline host 
parallel to the x axis, while the chain itself again lies along the z axis. The self-consistent 
equation for the dipole moments again takes the form ( [T2|) . but the diagonal elements of Q 
are given by ^ii(x, x') = d[diG{'s. — x'), where now 



^(x - x') 



A 1/2 

47re_Le|/ 



+ 

=11 



(22) 



For the case of interest, x = nc/z = 2;„z, x' = n'di = -2„/z, and one finds that Gxx{^ — x') = 

1 e^''^ 1 /111 /111 

QvvV^ ^ X ) ~ XT 1/2 1/2 u .13 and (x — x ) = ^75— 1/0 1/2 u .13 ■ The 



self-consistency condition determining the relation between uj and k can be written out, 
for all three polarizations, including only nearest neighbor dipole-dipole interactions, in the 
form 1 = — Stt^ i_r '4e ■ ^^^('^) cos/cd. Substituting in the values of Qu for the three cases, we 
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obtain 



2a^ Se^T ^^^'^ 



cos kd 



2a^ Seyy 1 



COS kd 



' = 35frrr-fci^<=-"- (23) 

Here Tj.j. = ry, Tyy = V^z = where Fy and are given by eqs. ffT7|) . Similarly, 
^Cxx = e(i^) — while 6eyy = 6ezz = e(i^) — e±- These equations can again be solved for 
k{uj) in the three cases, with or without a finite r, leading to dispersion relations with or 
without single-particle damping. 



C. Chain of Metallic Nanospheres in an External Magnetic Field 

Next, we turn to a chain of metallic nanoparticles in an external magnetic field, which 
we initially take to be parallel to the z axis. For such a system, we assume that the metal 
dielectric tensor is of the form ([T]), ([2]) and ([3]). The cases of a dilute suspension of metallic 



particles 21|, or of a random composite of ferromagnetic and non- ferromagnetic particles 24l |. 
have been treated previously. 

Once again, we take the chain of particles to lie along the z axis, with the n*'' particle 
centered at z„ = nd. The self-consistent equation for the dipole moments is still eq. f|T2|) . 
but now the elements of both Q and F are different from the case of an NLC host. For a 
chain of particles parallel to the z axis, Q is diagonal, with non-zero elements Qxx{^ — x') = 
Gyy{^ - x') = ^^i^^i^^^p , Gzz{y- - x') = - ^^j^^L^^^, |3 , whcrc Zn = ud aud we have assumed 
that the host has a dielectric constant equal to unity. The tensor F is also diagonal, with 
nonzero elements Ta = —1/3, i = 1, 2, 3. The quantity 6e = e{uj) — 1, where e{u) is now the 
dielectric tensor of the metallic particle. Using eq. ([1]) to evaluate this tensor, we obtain the 
following result for the tensor t = 6e[l — F^e]"^, to first order in the quantity A{uj), which 
is assumed to be small: 

txx tyy ^^^^(l ^ xx^^xx) 
txy tyx ^^xx^ xx^^xyi^^ ^ xx^^xx^ • (^^) 



Using these expressions, we can now write out the self-consistent hnear equations for the 
oscillating dipole moments and obtain dispersion relations for the modes. Once again, the 
longitudinal and transverse modes decouple. For the longitudinal modes, the self-consistent 
equation simplifies to 

^e{uJ) - 1 Pn' 



Pnz = 2a 



€{uj) + 2 ^ \z, 



Pn'z 

2-^ \y — yl \3' 



(25) 



This is the same as the equation for the longitudinal modes in the absence of a magnetic field 

and gives the same dispersion relation. For the transverse modes, the x and y components of 
the polarization are now coupled, and satisfy the equations [writing ^jj(x— x') = Gij{zn—Zn')] 



An 



Pnx Q ^ ^ [ixxGxxi^n ^n'^Pn' ,x ~l~ ixyGyyi^n Zn')P- 



47r 



Pny Q ^ ^ ^ {tyxGxxi^n ^n'^Pn' ,x ~^ ^yy^yyi^^n Zn')Pn',y 



(26) 



We can simplify these equations using the fact that Qxx{zn — ^n') = Syyi^n — Zn'), txx = tyy. 



and txy = tyx to obtain 



Pn,+ = - Y^^^-L - ^(^n - Zn')Pn',+ 



(27) 



and 



Pn, 



"^QiZn- Zn')Pn',-, (28) 

where t±^± — txx^itxy andp„,± = Pnx^iPny Thus, the equations for left- and right-circularly 
polarized waves are decoupled. 

To obtain explicit dispersion relations for left- and right-circularly polarized waves, we 
assume, as before, that p„ -t = p± exp{iknd — iut), and substitute the known forms for the 
quantities txx, txy, and Gxx{zn — Zn'), with the following result: 



1 = -2 



± 



?,A{u) 



E 



cos{nk±d) 



n=l 



(29) 



.e(a;) + 2 [e{u) + 2]\ 

In the special case where we include dipolar interactions only between nearest neighbors, 
this relation becomes 



1 = -2 



e{uj) - 1 ^ 3A{co) 



cos {k±d). 



(30) 



.e(c^) + 2 [e{uj) + 2]\ 
Since the frequency-dependence of both e(a;) and A{uj) is assumed known, these equations 
represent implicit relations between ou and k± for these transverse waves. 
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By solving for k±{uj) in eq. (129|) or f l30p . one finds tfiat left and right circularly polarized 
transverse waves propagating along the nanoparticle chain have slightly different wave vec- 
tors A;+ and for the same frequency u. Since a linearly polarized wave is composed of 
an equal fraction of right and left circularly polarized waves, this behavior corresponds to a 
rotation of the plane of polarization of a linearly polarized wave, as it propagates down the 
chain, and is analogous to the usual Faraday effect in a bulk dielectric. The angle of rotation 
per unit chain length may be written 

0M = ^[fc+M-A:_M]. (31) 

In the absence of damping, 6 is real. If r is finite, the electrons in each metal particle will 
experience damping within each particle, leading to an exponential decay of the plasmons 
propagating along the chain. This damping is automatically included in the above formalism, 
and can be seen most easily if only nearest neighbor coupling is included. The quantity 

e{u) = 9i{u) + i92{u;) (32) 

is then the complex angle of rotation per unit length of a linearly polarized wave propagating 
along the chain of metal particles. By analogy with the interpretation of a complex ^ in a 
homogeneous bulk material, Re6'(ci;) represents the angle of rotation of a linearly polarized 
wave (per unit length of chain), and lm6{u}) as the corresponding Faraday ellipticity - i. e., 
the degree to which the initially linearly polarized wave becomes elliptically polarized as it 
propagates along the chain. 

For a magnetic field perpendicular to the chain (let us say, along the x axis), the elements 
of the matrix t become 

txx = tyy = t^^ = 3[e{u) - l]/[e(w) + 2] 

tyz = -tzy = -3tA{w)[e{uj) - l]/[e{io) + 2f , (33) 

with other elements equal to zero. The transverse waves polarized parallel to the x axis are 
unaffected by the magnetic field, and are described by the equations 

Pnx 3"*^ ^xx ^ ^ Qxxi^^n ^n'^Pn' ,x- (^^) 
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The y and z polarized waves, with components pny and pm-, are coupled, however, and satisfy 



Pny 



Pn 



2 ^ ^ {tyyQyyi^Zn Zn'^Pn' ,y ~\~ tyzQ zzi^Z^ Z^i^Pn' ,z\ 

^ ^ [tzyGyyi^n Zn')Pn',y + ^zzOzzi^n ^nOPra'jz] (^5) 

The tensor ^ is still diagonal, with the same nonzero elements as in the case of B parallel 
to the chain. Assuming propagating waves of the form pny = poy e^p{inkd — iut), pnz = 
Po2exp(mA;(i — iut), we find that the amplitudes poy and poz satisfy the equations 



POy 



POz 



^ ^ [^yySyyi ^n' )poy + tyzQzz{-Zn')poz] exp{iknd) 



— ^ [-tyzQ: 



-Zn')poy + tzzQzz{-Zn')poz] exp{ikn'd). 



(36) 



In the special case where only nearest neighbor interactions are included, these equations 
simplify to 



POz 



Sir 

[tyyGyyid)poy + tyzQzz{d)pQz\ COs{kd) 

Stt 

[-tyzQyy{d)poy + tzzG zz{d)poz] COs(/crf). 



(37) 



Substituting in the explicit forms of Qyy and Qzz, we find that these equations take the form 



Poy 



POz 



a 
rf3 



a 



^tyyPQy + ^yzPOZ 

'2 4 

-^tyzPoy + -^tzzPOz 



COS kd 

COS kd. 



(38) 



If we solve the pair of equations ( 136|) or ( 138|) for poy and poz for a given value of k, we 
obtain a nonzero solutions only if the determinant of the matrix of coefficients vanishes. For 
a given real frequency u, there will, in general, be two solutions for k{u) which decay in 
the +z direction. These correspond to two branches of propagating plasmon (or plasmon 
polariton) waves, with dispersion relations which we may write ki{u) and k2{u). The fre- 
quency dependence appears because both tyz and tyy depend on u [through e(co') and A(w)]. 
The corresponding solutions {poy, poz) are no longer linearly polarized but will instead be 
elliptically polarized. However, unlike the case where the magnetic field is parallel to the 
z axis, the waves are not circularly polarized, and the two solutions are non-degenerate. 
Because A{u) is usually small, the ellipse has a high eccentricity, and the change in propa- 
gation characteristics due to the magnetic field will usually be small for this magnetic field 
direction. 
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III. NUMERICAL ILLUSTRATIONS 



As a first numerical example, we calculate the plasmon dispersion relations for a chain of 
spherical Drude metal particles immersed in an NLC We consider two cases: liquid crystal 
director parallel and perpendicular to the chain axis, which we take as the z axis. For ey 
and e i , we take the values used in Ref. 25|. These are taken from experiments described in 
Ref. j26|, which were carried out on the NLC known as E7. For comparison, we also show 
the corresponding dispersion relations for an isotropic host of dielectric constant which is 
arbitrarily taken as |e|| + |e_L = 2.5611. The results of these calculations are shown in Figs. 
1 and 2 in the absence of damping (r — t- cxd in the Drude expression). As can be seen, both 
the L and T dispersion relations are significantly altered when the host is a nematic liquid 
crystal rather than an isotropic dielectric; in particular, the widths of the L and T bands 
are changed. When the director is perpendicular to the chain axis, the two T branches are 
split when the host is an NLC, whereas they are degenerate for an isotropic host, or an NLC 
host with director parallel to the chain. 

Next, we turn to the effects of an applied magnetic field on these dispersion relations for 
the case B||z. The L waves are unaffected by a magnetic field, but the T waves are split 
into left- and right-circularly polarized waves. To illustrate the predictions of our simple 
expressions, we again take ajd = 1/3, and we assume a magnetic field such that the ratio 
Uc/ujp = 3.5 X 10~^. For a typical metallic plasma frequency of ~ 10^^ sec~^, this ratio 
would correspond to a magnetic induction _B ~ 2 T. We consider both the undamped case 
(r — 7- oo) and the damped case [upT = 100). Using these parameters, the dispersion relations 
for the two circular polarizations are shown in Fig. 3 both with and without single-particle 
damping. The sphtting between the two circularly polarized T waves is not visible on the 
scale of the figure. 

We also plot the corresponding rotation angles 6d for a distance equal to one interparticle 
spacing in Fig. 4. When there is no damping, 6 diverges near the edge of the T bands, but 
this divergence disappears for a finite r (e. g. UpT = 100, as shown in in Fig. 4). In this 
case, Re9{uj)d never exceeds about 0.001 rad per interparticle separation, corresponding to 
a rotation of less than 0.1° over this distance. Im6{u)d is also small, showing that a linear 
incident wave acquires little ellipticity over such distances. Over a distance of thirty or 
so interparticle separations, a linearly polarized transverse wave would typically rotate by 
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only 1-3°. Since theory and experiment suggests that the wave intensity typically has an 
exponential decay length of no more than around 20 interparticle spacings |27|, the likely 
Faraday rotation of such a wave in practice will probably not exceed a degree or two, at 
most, even for a field as large as 2 T. Thus, while the rotation found here is likely to be 
measurable, it may not be large, at least for this simple chain geometry with one particle 
per unit cell. The present expressions also indicate that 6 is very nearly linear in B, so a 
larger rotation could be attained by increasing B. 

As can be seen, 9 depends strongly on frequency in the case of zero damping, but less 
so at finite damping. In the (unrealistic) no-damping case, at the very top and the very 
bottom of the plasmonic band, only one of the two circularly polarized waves can propagate 
down the chain. Because this filtering occurs only over a very narrow frequency range (of 
order Uc/cOp), and because this calculation assumes no damping, it would be quite difficult 
to detect a region where only one of the two circularly polarized waves can propagate. 

Finally, we mention the case where B ± z for the same parameters as in the parallel 
case. The effect of a finite B causes two non-degenerate dispersion relations (one an L and 
the other a T wave) for the perpendicular case to become mixed. We have not computed 
the rotation angles for this perpendicular case, but we expect that they would be similar in 
magnitude to the values shown in Fig. 4. 

IV. DISCUSSION 

The calculations and formalism presented in the previous section leave out several effects 
which may be at least quantitatively important. First, in our numerical calculations, but 
not in the formalism, we have included only nearest neighbor dipolar coupling. Inclusion of 
further neighbors will quantitatively alter the dispersion relations in all cases considered, as 
well as the Faraday rotation angle when there is an applied magnetic field, but these effects 
should not be very large, as is already suggested by the early calculations in Ref. [lO| for 
an isotropic host. Another possible effect will appear when a/d is significantly greater than 



1/3, namely, the emergence of quadrupolar and higher quasistatic bands 12|. These will mix 
with the dipolar band and alter its shape. For the separations we consider, this multipolar 
effect should be small. Also, even if a > d/3, the plasmon dispersion relations will still be 
altered by an NLC host or by an applied magnetic field in the manner described here. 
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The present treatment also omits radiative damping. In the absence of a magnetic field, 
such damping is known to be small but non-zero in the long- wavelength limit, but it becomes 
more significant when the particle radius is a substantial fraction of a wavelength. Even at 
long wavelengths, radiative damping can be very important at certain characteristic values 
of the wave vector [l^. We have not, as yet, extended the present approach to include such 
radiative effects. We expect that, just as for an isotropic host in the absence of a magnetic 
field, radiative effects will further damp the propagating plasmons in the geometries we 
consider, but will not qualitatively change the effects we have described. 

For the case when the host is an NLC, the present work oversimplifies the treatment of 
the NLC host by assuming that the director field is uniform, i. e., position-independent. In 
reality, the director is almost certain to be modified close to the metal nanoparticle surface, 
i. e., to become nonuniform, as has been pointed out by many authors^ 



The effects of 



such complications on the optical properties of a single metallic particle immersed in an 
NLC have been treated, for example, in Ref. and similar approaches might be possible 
for the present problem also. 

To summarize, we have shown that the dispersion relations for plasmonic waves prop- 
agating along a chain of closely spaced nanoparticles of Drude metal are strongly affected 
by external perturbations. First, if the host is a uniaxially anisotropic dielectric (such as 
a nematic liquid crystal), the dispersion relations of both L and T modes are significantly 
modified, compared to those of an isotropic host, and if the director axis of the NLC is 
perpendicular to the chain, the two degenerate transverse branches are split. Secondly, if 
the chain is subjected to an applied magnetic field parallel to the axis, the T waves undergo 
a small but measurable Faraday rotation, and also acquire a slight ellipticity. A similar ellip- 
ticity develops if the magnetic field is perpendicular to the chain, but its effect will likely be 
more difficult to observe, because the magnetic field couples two non-degenerate branches 
of the dispersion relation. All these effects show that the propagation of such plasmonic 
waves can be tuned, by either a liquid crystalline host or a magnetic field, so as to change 
the frequency band where wave propagation can occur, or the polarization of these waves. 
This control may be valuable in developing devices using plasmonic waves in future optical 
circuit design. 
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FIG. 1: Calculated dispersion relations uj{k) for plasmon waves along a chain of metallic nanopar- 
ticles, in the presence of an NLC host. We plot oj/ujp, where ujp is the plasma frequency, as a 
function of kd, where d is the distance between sphere centers. Green and blue (x's and +'s): L 
and T modes for a chain embedded in an NLC with director parallel to the chain. The NLC is 
assumed to have principal dielectric tensor elements ey = 3.0625 and = 2.3104 parallel and 
perpendicular to the director, corresponding to the material known as E7. In this and all subse- 
quent plots a/d = 1/3, where a is the metallic sphere radius. Also shown are the corresponding L 
and T dispersion relations (black and red solid lines, respectively) when the host is isotropic with 
dielectric constant e/j = 2.5611 = iey + |e_L. 
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FIG. 2: Same as Fig. 1 except that the director of the NLC is perpendicular to the chain of metal 
nanoparticles. The frequencies of the L modes (asterisks, in green) and T modes (+'s and x's, 
shown in dark and light blue), divided by the plasma frequency Wp, are plotted versus kd. The 
NLC has the same dielectric tensor elements as in Fig. 1. Also shown are the corresponding L 
(solid black) and T (solid red) branches for an isotropic host with e/j = 2.5611. Note that the T 
branches which were degenerate in Fig. 1 are split into two branches in this NLC geometry. 
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FIG. 3: Blue symbols (+ and x symbols): the dispersion relations for left and right circularly polar- 
ized T plasmon waves propagating along a chain of metallic nanoparticles, assuming no damping 
and a magnetic field parallel to the chain. We assume also that coc/^p = 3.5 x 10"^, where 
Uc = eB/mc is the cyclotron frequency. Green + and x symbols: same dispersion relations, but 
assuming single-particle damping corresponding to cOpT = 100. In both cases, the splitting be- 
tween left and right circularly polarized waves is not visible on the scale of the figure, so that for 
both undamped and damped dispersion relations there appear to be single sets of asterisks. For 
Up = I. X 10^^ sec"-*^, the chosen ujc/tdp corresponds to about 2 T. 
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FIG. 4: Solid curve: Faraday rotation angle 9d per interparticle spacing as a function of frequency, 
assuming r — )• oo, and oJc/^p = 3.5 x 10^^. Blue and red (+ and x) symbols: Real and imaginary 
parts of the Faraday rotation angle per interparticle spacing, 9d, assuming UpT = 100 and the ratio 
b^d^v ~ "^-^ ^ 10^^. Note that there are two different vertical scales for undamped and damped 
rotation angles. 
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